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Abstract 

We study a symmetry broken harmonic oscillator in contact with a heat bath characterized by a 
fixed temperature. The overdampped system is solved exactly to show symmetry broken directed 
transport raising the question whether fluctuation-dissipation relation does guarantee equilibrium. 
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Consider two bodies of same mass but of different damping (due to difference in shape size 
etc) are connected by a spring and the system is in contact with a heat bath of temperature 
T. The equations of motion are the following 
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In the above system the coordinates x\ and x 2 are the center of mass of the bodies 1 and 2 
where x\ > x 2 and the hard core volume exclusion would be applied (later) by taking the 
limit V — > infinity. The spring constant is a, < rji >= and < r)i(t)r)j(s) >= 5ij5(t — s) 
ensures the system to be at a temperature T according to the standard fluctuation dissipation 
relationship where k B is the Boltzmann constant. 

Taking the overdamped limit and moving to the internal coordinate Z = x\ — x 2 and 
centre of mass X = X1 ~ X2 we can re write the above system as 

• a(2-/3) t V(2-/3)dQ(-Z) 

Here < iz(t)i z (s) >= 2Tk B ^5(t - s) and < &(t)&(s) >= Tk B ^^5(t - s) 
The distribution of Z follows from Eq.3 as 
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The hardcore repulsion term at the limit V tending to infinity will suppress the distribution 
on the negative Z space and that will result in a nonzero < Z >. The average velocity of 
the cm being < X >= — 2 ^-p) < Z > + 2 ^f^ J* e °° S(Z)P(z)dz will now be nonzero where e 
is an infinitesimal positive constant. The second term in the expression for < X > is zero 
for a nonzero e where e = is not accessible because of excluded volume i.e. the Heaviside 
step function Q(—Z) or equivalently by the vanishing of the P(Z) at Z < 0. 

Let us note that no contribution from the second term in the expression of < X > may be 
understood from a mechanistic view point. All the interactions between two bodies under 
consideration such as hard core collision or spring interaction are themselves momentum 
preserving. The resultant motion of the CM is due to different damping experienced by two 
bodies on the flight due to the interaction with the environment. Thus, the requirement 



according to the Newton's first law is fulfilled where there is an external force (damping) that 
causes nonzero velocity of the CM of the system. Now, the hard core collision (infinite V) is 
ideally a momentary event when the velocity of the colliding particles are zero. Just before 
and after the collision the particles have velocities. So, the hard core collision should remain 
a momentum preserving event and as a consequence there should not be any contribution 
to the velocity of the CM from the hard core collision which justifies the vanishing of the 
contribution to < X > from the hard core collision term. Another point to note that, in 
the above calculations, the qualitative results are practically independent of the strength of 
the thermal noise (rji) so long as rji and r]2 are uncorrelated (which practically led to the 
presented calculations from a previously otherwise considered model To bring such a 
system to equilibrium would there be a correlation established in t^s? If so, it goes somewhat 
against the existing concept of representation of a heat bath in equilibrium. 

This is something extraordinary! Extraordinary because the system is not subjected to 
any external force simply the symmetry has been broken by the difference in the damping 
constants and the strength of the Gaussian noise being proportional to the damping con- 
stants keeps the temperature seen by the two bodies making up the system constant. In 
the standard form of a Langevin equation there always is the scope to keep an external 
force felt by the particle and a system can always move in a constant temperature heat bath 
under the action of this external force showing some fluctuations on its deterministic path 
of motion. But here, there is an internal field felt by the system and still it moves being 
in contact with a heat bath characterized by a constant temperature as per the demand of 
the fluctuation-dissipation relation. But, in equilibrium a system cannot move direction- 
ally which is in conflict with the concept of the stationary probability distribution. So, 
the question naturally arises as to what extent the derivation of the fluctuation-dissipation 
relation for a single particle in equilibrium with a heat bath applies for a system of particles 
having an internal time scale (1/a). In the derivation of the FDT its generally taken into 
account that the heat bath equilibrates very quickly [2] and that too without any reference 
to a particular time scale with respect to which it should be quickly. But when, there is a 
particular characteristic time scale involved should it be defined differently if we assert that 
there is nothing wrong with the overdamped limit? 
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I. ZWANZIG MODEL 



Let us derive the model Eq.l and Eq.2 following the method of Zwanzig's to show that 
under certain conditions where the couplings of the two bodies with the environmental 
degrees of freedom are uncorrelated the cross terms (in momentum) in Eq.l and 2 can be 
avoided and Eq.l and 2 retain their shape and do not become a standard Ornstein-Uhlenbeck 
one. Consider the Hamiltonian of the system as 

I \ it/ j 

In the above expression, g^s are the bath degrees of freedom. The dynamics that results 
from the above mentioned equation are 

dV(X 1 ,X 2 ) , ^ / 7h*i72^2 
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and 



Pi + ufqi = 7 H Xi + 7 2i X 2 (9) 
Following the standard procedure Eq.7 and 8 can be rewritten as 
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In the above equations fi(t) and f 2 (t) are initial condition dependent. Since there is a sum 
over the bath degrees of freedom in the cross terms containing ju^i, under the condition 
that the 7^ and 7^ are uncorrelated this cross term should vanish making the equations 10 
and 11 take the shape of Eq.l and 2. So, the form of Eq.l and 2 makes sense where we 
have not taken into account the cross terms as is done in a standard Ornstein-Uhlenbeck 
process. Now, the results shown above are a consequence of the structure of Eq.l and 2 and 
that holds at least under the situation when the coupling of the two bodies (connected by 
spring) to the environmental degrees of freedom is uncorrelated. 

Zwanzig model although is a widely accepted microscopic picture leading to the Langevin 
dynamics. However, its physically far from reality because the long range strong coupling 



between all the degrees of freedoms. Nevertheless, it at least, when put in the context of the 
introductory section's results, indicates that the uncorrelated coupling of the system with 
the degrees of freedom of the bath is something leading to the breakdown of second law here 
and breakdown of second law is something totally unacceptable. So, the indication made 
in the introductory section that to restore equilibrium i.e. restore second law a correlated 
forcing of the system would be necessary is at least getting some support following this 
Zwanzig-model calculations. 

If we think physically, how would the equilibrium be established? Consideration of inher- 
ent uniformity (or symmetry) of an equilibrium state when put together with the symmetry 
broken dimer is leading to all these problems. So, if the dimer is symmetry broken (shape 
and size of the two bodies are different) then their local atmosphere will also get symme- 
try broken by having variations in density of particles etc. This is a situation where the 
change in the local environment of the two bodies making the dimer would restore the bro- 
ken symmetry making the damping of two very different objects the same. Now, this whole 
locally varied environment of the dimer has to come to terms with the global uniform (as 
equilibrium normally demands) state and that remains an issue. The possibility of a global 
symmetry breaking cannot be ruled out even in equilibrium and probing its dynamics since 
the insertion of the dimer into an uniform equilibrium system might be interesting. 



II. NON-SINGULAR COLLISION POTENTIAL 



It may appear that the result shown above is a consequence of a hard core potential 
only and to counter that in this section we would consider a nonsingular potential for the 
collision of the two bodies constituting the dimer {3]. Instead of considering the ©(— Z) 
a Heaviside step function, let us take the potential as Q(Z) = —a(z — zo)forZ < Zq and 
Q(Z) = OforZ > Zq. The discontinuity at the point Zq is essential as the constituents 
of the dimer are just coming in contact with each other at this gap. Moreover, we would 
further consider that the gradient a of this repulsive potential is steep enough to ensure non- 
passage of the classical particles past each other to make Z negative. So, a and are adjustable 
parameters at our hand. With such a potential the dynamics on internal coordinate and the 
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CM now takes the form 



Z = - 



d_ 

dZ 



a(2-(3) z2 _2-(3 



2(1-/3) 
X - 



1-/3 



a(Z — Z )f orZ <z 



a/3 



-Z + 



2(1-/3) 2(1-/3) 
The probability distribution for the Z would be 



0>forZ<Z 



r .... , aZ 2 a(Z — Z ) 

P(Z) = exp ( -— - + 1 0) 

ZKb-L KbJ- forZ<Z 



Now, we can write the average velocity of the CM as 



<X >=- < 



a/3Z 
2~p 



forZ<Z 



> + < 
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forZ<Z 
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(12) 
(13) 

(14) 

(15) 



The averages on the right hand side are now for the Gaussian probability distribution 



exp 



2k B T 



only. Definitely these two terms on the right hand side of the above equa- 
tion would not cancel in general being the averages of two very different quantities with a 
lots of free parameters involved. Note that, by the consideration of a being large enough we 
are doing away with any probability of Z < and the integration limit is practically from 
to oo while averaging. Thus, any objection considering the delta function singularity (or 
equivalently infinitely discontinuous potential) only resulting in the CM velocity does not 
stand. So long as the excluded volume is there the result is there. 

Let us try to understand why the resultant velocity of the CM is practically independent 
of the type of the repulsive potential. The repulsive excluded volume potential and the 
spring potential can actually be two very different things giving the scope of having many 
tunable free parameters to do away with any cancelation of the positive and negative parts of 
the CM velocity resulting in zero resultant velocity. Physically, it means that you can always 
control the range of the flight of the two bodies making the dimer before they collide and how 
they collide and the external forcing resulting in the CM motion practically works during 
the flight. Thus, issue whether the fluctuation-dissipation relation guaranties equilibrium, 
now, more certainly does indicate that it does not always do that. 



III. HOW TO EXTRACT ENERGY 



In order to design an engine to extract energy from a quasi one dimensional bath one 
has to simply introduce sush an asymmetric dimer into the system. As soon as the dimer is 
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introduced the system turns nonequilibrium and starts its journey to the next equilibrium 
state where the density of the environment would be different on the two sides of the dimer. 
Actually, one can imagine that the dimer moves to achieve this density varied state by 
changing the effective volume on either sides of it. So, long as the required density veriation 
is achieved the dimer keeps moving to adjust the volume. Note that, here we are considering 
the quasi ID heat bath is in contact with a bigger heat bath (atmosphere say) to keep its 
temperature constant. Now, once the dimer has stopped following reverse density gradient 
being established which is strong enough to stop it, one simply has to remove the dimer and 
allow the system to come back to the initial state. If the energy cost in the introduction 
and removal of the dimer is less than the energy extracted from its motion there is always a 
gain in energy. Most important to note that here one is extracting the energy directly from 
the bath unlike in Brownian ratchet models where work is done by the external forcing. 
So, the proposed mechanism is a completely novel one and is in contrast with the standard 
paradigm of Brownian Ratchetting. 



IV. ORNSTEIN-UHLENBECK STRUCTURE 



The model we are considering here in 2D differs from the structure of a standard Ornstein- 
Uhlenbeck (OU) form in the sense that it does not include the cross terms of velocities in the 
eq.l and eq.2. A possible justification for that has been given in the section Zwanzig-model. 
Here we will note that there are certain problems with the OU form of this model as well 
because of the singular nature of the matrix of the elastic force. The standard OU form for 
our model is 

(-* a )( X ') + hu (16) 
D 2 I j \x 2 J ya -a J \x 2 J \r) 2 J 

where we have considered hardcore collision potential of practically no width (i.e delta forcing 
at collision). This can be taken into the form of a standard OU process as 

2 

-V, • », ; .Y 'ir, (17) 

i 

where 

' I " II I (18) 



One can readily find out the eigenvalues of the matrix T which are 



A 2 





q(2 + Pi + D 2 - P) 
I -f3- D 1 D 2 



(19) 
(20) 



The zero eigenvalue is definitely a problem here because for the system to reach a stationary 
state the eigenvalues have to be greater than zero. The surprising thing is that in this whole 
calculation we are not considering a strong external drive which basically takes the system 
out of the linear response regime where the fluctuation-dissipation relation does not in any 
way hold. Any strong interaction the system encounters apart from that with the bath 
degrees of freedom is through the collision of the two bodies to ensure volume exclusion and 
that is something pretty much naturally occuring everywhere. If this collision is taking the 
system out of the linear response regime then it shoudl always happen and the credibility 
of the fluctuation-dissipation relation in its standard form is questionable. 
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